Factorization Properties in the 3D Edwards-Anderson Model 
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Starting from the study of a linear combination of multi-overlaps which can be rigorously shown 
to vanish for large systems we numerically analyze the factorization properties of the link-overlaps 
multi-distribution for the 3D Gaussian Edward- Anderson spin-glass model. We find evidence of a 
pure factorization law for the multi-correlation functions. For instance the quantity 

< <3i2 > - < Q12Q34 > 



< 



tends to zero at increasing volumes. We also perform the same analysis for the standard overlap for 
which instead the lack of factorization persists increasing the size of the system. The necessity of a 
better understanding of the mutual relation between the two overlaps is pointed out. 



PACS numbers: 05.50.+q, 75.50.Lk 



The structure of the low temperature phase for the fi- 
nite dimensional spin glass is among the most interesting 
and yet unsettled problems in condensed matter. After 
more than thirty years its main issues remain unsolved 
not only from the mathematically rigorous point of view 
but also within the theoretical physics perspective. In 
particular it is not clear what is the quenched probability 
distribution of the spin-overlap for large systems and low 
temperatures and different pictures have been proposed: 
the Replica Symmetry Breaking [IJ which describes the 
distribution similarly to the mean field one, the Droplet 
picture |2j which claims that the support of the distribu- 
tion shrinks to a point (becomes trivial) when the system 
size grows to infinity (see also 0. [J for other possible 
pictures and different perspectives). 

In this paper we investigate the structure of the multi- 
distribution of the overlaps among replicated samples of 
the Edwards- Anderson (EA) model with Gaussian near- 
est neighbor couplings. The sp in g lass quenched mea- 
sure is indeed described (see [y, [Tl| for instance) by an 
infinite family of probability distributions which repre- 
sent its equilibrium state: the distribution of the single 
overlap P{Q\,2) related to the internal energy of the sys- 
tem but also those involving more than two replicas like 
P(Qi,2, Qi,s)) P{Q\,2, 63,4) related to the specific heat 
etc. Since in the quenched measure the different copies 
are taken with the same frozen disorder the random vari- 
ables <2i, m are not independent and their joint distribu- 
tion does not factorize on products of the single one at 
finite volumes: P(Qi, 2 , 61,3) # ^(61,2)^(61,3). 

With this work we address precisely the question if 



such a factorization may occur when the thermodynamic 
limit is reached and we find strong numerical evidence for 
a positive answer in the link overlap case and a negative 
one for the standard overlap. 

Previous studies on the model 0, || have concentrated 
their attention on the single overlap distribution and 
claims were made about its triviality (for different per- 
spectives see also 0, 0, 0). We independently re- 
produced those numerical results and extended them to 
larger sizes (up to L = 12). In our opinion they can- 
not distinguish among the trivial or non-trivial picture. 
In particular the data for the link overlap do not rule 
out a limiting distribution with two peaks and a plateau 
among them. That fact together with a complete lack of 
rigorous results make the factorization properties of the 
multi-overlaps a completely open matter and motivate 
our investigation. 

Our departing point are the factorization rules found 
for the Edwards- Anderson model in terms of its link- 
overlap. In a box of side L the link overlap between two 
spin configuration a and r is defined as 



far) 



(1) 



where the sum runs over all couples of sites which are 
connected by a random bond (usually the nearest neigh- 
bor sites). In |(| and Q it was rigorously proved that 
when the size of the system grows to infinity the linear 
combination 



<Ql,2 



4Ql,2<22,3 + 3Qi, 2 Q 3 ,4 > 



(2) 
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is vanishing except possibly on isolated temperatures 
where phase transitions may occur. Here the brackets 
denote the quenched measure (i.e. the successive compu- 
tation of the expectation w.r.t. the Boltznmann-Gibbs 
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2032 
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TABLE I: Parameters of the simulations: system size, number 
of sweeps used for thermalization, number of sweeps during 
which observables were measured, number of disorder sam- 
ples, number of f3 values allowed in the parallel tempering 
procedure, temperature increment, minimum and maximum 
temperature values. 

measure on replicated samples followed by the average 
over the Gaussian disorder) and the subscripts on the 
Q's label the different real replicas. 

The same relation Q was originally found to hold for 
the square of the standard overlap 




in the SK model within the RSB pict ure and related to a 
property of replica equivalence [j^llOf . See also 0, El 0] 
for its rigorous derivations. 

Although the vanishing of © emerged first in the mean 
field picture its validity alone cannot distinguish between 
the different scenarios proposed (see the discussion in |(| ) . 
In particular it cannot distinguish between the peculiar 
ultrametric factorization rule like the one proposed for 
the SK model or a pure factorization rule of the multi- 
overlap distributions. 

In order to test the factorization properties we chose to 
consider generic linear combinations of the above mono- 
mials with coefficients whose sum is zero: 

g(a) =<Q?, 2 -"Qi,2Q2,3 + (a-l)Qi,2Q3,4> ■ (4) 

From we know that the former expression is close to 
zero in a = 4 with finite volume correction of size L~ 3 . 
Away from a = 4 a link-overlap multi-distribution with 
a pure factorization law between replicas would predict a 
progressive squeezing to zero (at increasing volumes) of 
the line g(a) for all a . Instead a non-factorizing link- 
overlaps distribution would be compatible with the per- 
sistence away from zero of the line angular coefficient 

m = < Ql,2<33,4 - Ql,2<92,3 > (5) 

and intercept 

n =< Q\ 2 - Qi, 2 <2 3 ,4 > ■ (6) 

We performed the analysis for both the standard over- 
lap and for the link overlap. Before going to the detailed 



descriptions of the model and the illustration of the re- 
sults it is worth to mention the different roles played by 
the two. The standard overlap is historically the first 
proposed observable in the study of the spin glass phase 
being directly related to the original Edwards- Anderson 
order parameter. Its widespread use is moreover due also 
to the fact that in the SK model its distribution carries 
the whole information of the thermodynamic properties. 
Undoubtedly it is a very interesting quantity to be stud- 
ied for the low temperature phase of general disordered 
spin systems; nevertheless we want to point out that each 
spin glass model has its own natural observable which, 
in a Gaussian model, is given by the covariance of its 
Hamiltonian. An easy computation |(| shows that in the 
EA model such a covariance is the link-overlap in terms 
of which all the thermal observable can be expressed like 
internal energy, specific heat etc, and rigorous results can 
be established like stochastic stability Q, Q ; see also [l^ 
for its use in rigorous results from a different perspec- 
tive. It is also interesting to observe that the standard 
overlap between two spin configuration changes propor- 
tionally to the volume of the different spins in the two 
configurations; instead the link overlap changes like the 
surface. 

Recalling the basic definitions: we consider the Gaus- 
sian Edward- Anderson model , defined by the Hamil- 
tonian 

H(J,a) = -^2 Ji^+e^iO-i+e^ (7) 

i£A fj.=x,y,z 

where i is a site of a 3-dimensional cubic lattice A (|A| = 
L 3 ), e M is the versor in the /i direction (with /x = x,y, z), 
Cj = ±1 are Ising spin variables and Ji t i+ e , are Gaussian 
random variables with zero average and unit variance. 

We performed numerical simulations by using the Par- 
allel Tempering (PT) algorithm to facilitate equilibra- 
tion. We used periodic boundary conditions and investi- 
gated lattice sizes up to L = 12. For every size we sim- 
ulated at least 2032 realizations of the couplings. Other 
parameters of the simulations are reported in Table {IJ. 
The allowed temperature range (assuming T c ~ 1.) was 
approximately 0.5T C < T < 2.3T C and we used up to 
37 temperatures in the PT procedure. We tested ther- 
malization by checking the symmetry of the probability 
distribution for the standard overlap under the transfor- 
mation q — ► —q. 

If not otherwise stated in the sequel we always plot 
the same quantities relative to the two overlap with the 
same scale on the y-axis, in order to let better appreciate 
analogies and differences between the two. 

Fig. (JTJ shows the plot of relation in Eq. © as a func- 
tion of the temperature; we see on the left side that the 
vanishing of the linear combination is well reproduced 
by numerical simulation for the link overlap (even for 
the small size L — 4 the finite size correction is less than 
5 ■ 10~ 3 ). It is interesting to see that also for the squared 
standard overlap (right side) the relation is satisfied (in 
agreement with ll3j), even if there are no rigorous argu- 
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ment which support it. We notice moreover that finite 
size corrections are larger than those for the link overlap. 

In Fig. J2J we show the different lines g(a) as the 
volume increases. For all the sizes the relation is close 
to zero for a = 4. When the system size increases we 
observe that the lines tend to flatten, which is a possible 
signal of triviality. We may observe that the flattening is 
much more evident for the link overlap compared to the 
standard overlap. 

To investigate further the factorization and triviality 
matter, we measure the angular coefficient (0 and the 
intercept © and see how they scale w.r.t. to the system 
size L. We find that both have some tendency to decrease 
and we fit the data with a law of the type y = a + c ■ L 5 
for different values of a and measuring the relative chi- 
square. We find that the minimum chi-square for the an- 
gular coefficient and intercept relative to the link-overlap 
is reached for a = which supports the factorizing pic- 
ture. More precisely the data gives a normalized x 2 = 1-5 
for a = and the \ 2 value is increasing with a, being al- 
ready of O(10) for a = 0.001. On the other hand the 
same analysis for the standard overlap showed that the 
chi-square is basically constant for all values of a close 
to zero. It has a normalized \ 2 — for all values 

in the interval [0,0.001]. These results suggests that the 
data we found for the standard overlap do not allow to 
distinguish among a factorizing picture (a — in the 
thermodynamic limit) and a non-factorizing (o ^ 0). 

Finally we analyzed the normalized angular coefficient 
m* = mj < Q\ 2 > and intercept n* = n/ < Q\ 2 > 
which are adimensional quantities normalized with their 
typical values (they are slightly different for the two over- 
laps). The data reported in log- log scale in Fig.© show 
a clear factorization tendency for the link overlap multi- 
distribution while the standard overlap seem to lack the 
same property. The chi-square analysis confirms all that. 
This works shows that, within the current reachable lat- 



tice sizes, the multi distribution for the link overlap obeys 
to a pure factorization law while the standard overlap 
lacks the same property. A similar factorization has a 
very clear meaning: the random variables Qi t , n become 
independent in the infinite volume limit with respect to 
the quenched measure. This fact is certainly a form of 
triviality and is observed for instance in the Curie- Weiss 
model where the spin variables <Ji&j become independent 
in the infinite volume limit with respect to the Boltzmann 
measure (see ||). We want to stress nevertheless that the 
factorization triviality that we found is compatible with 
both the usual trivial or nontrivial pictures because it 
doesn't imply anything on the single overlap distribution 
but only on the multi overlap ones. 

We observe, finally, that the different behavior of the 
two overlaps requires a better understanding of their rel- 
ative fluctuations: does the standard overlap distribution 
concentrate at fixed link overlap? And if it does is the 
functional relation among the two a one-to-one function? 
The two questions can be trivially answered in the SK 
model where the link overlap is the square power of the 
standard overlap. In a finite dimensional model like EA 
the relation among the two is much more involved and the 
problem of their relative fluctuation has been numerically 
addressed in at T = and more recently in 0| for 
a model with soft constraints. Their results point to the 
direction of a vanishing relative fluctuation. The results 
we found on this paper would be compatible either with 
a non vanishing relative fluctuation or with a multivalued 
functional relation among the two random variables that 
could account for the difference in factorization behavior. 
We plan to investigate further those two open questions 
both numerically and analytically. [2(J 
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results are shown in Fig.JIJ. What we observe is the re- 
markable coincidence of the np and mp for the two dif- 
ferent overlaps. Introducing the correlation coefficient p 
and observing that np = I — p the identical behavior of 
p for the two overlaps suggests a high mutual correlation 
among them. The fact that p stays basically constant in 
the observed volume range does not give much informa- 



tion on the factorization of the multi overlap because it 
may happen that the variance at the denominator shrinks 
to zero like the covariance at the numerator. Instead the 
coefficient k normalized with the second order moment 
shows the same behavior of m* and n* and confirms our 
conclusions. 
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FIG. 1: Plot of relation as a function of the temperature for different system sizes (see the legend). On the left data for 
the link overlap, on the right data for the standard overlap. 
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FIG. 2: Plot of the function g(a) versus a for the fixed temperature T — 0.6: link over lap (left), standard overlap (right). 
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FIG. 3: The absolute value of the normalized angular coeffi- 
cient m* (black circles for the link overlap, green diamond for 
standard overlap) and intercept n* (red squares for the link 
overlap, blue triangles for the standard overlap) versus the 
lattice size L. The temperature is T = 0.6. Lines are guide 
to the eyes. 
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FIG. 4: A different normalization of the data (see the legend). 
Here x = Q12, y = Q34 (left panel) and x — g 12 , y — 934, 
(right panel), V^(a;) = (x 2 ) — (x) 2 — V(y) denotes the variance 
and C(x,y) — (xy) — (x)(y) denotes the covariance. 



